Entangling a nanomechanical resonator and a superconducting microwave cavity 
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We propose a scheme able to entangle at the steady state a nanomechanical resonator with 
a microwave cavity mode of a driven superconducting coplanar waveguide. The nanomechanical 
resonator is capacitively coupled with the central conductor of the waveguide and stationary entan- 
glement is achievable up to temperatures of tens of milliKelvin. 
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I. INTRODUCTION 

Entanglement is one of a number of inherently quan- 
tum phenomena that it is hoped will soon be observable 
in macroscopic mechanical systems Aside from the 
interest in studying quantum mechanics in a new regime, 
entanglement may be used as part of read-out schemes 
in quantum information processing applications. Meth- 
ods for entangling a nanomechanical resonator with a 
Cooper pair box [2[, or an optical mode [3], for entan- 
ling two charge qubits or two Josephson junctions 
via nanomechanical resonators, and for entangling two 
nanomechanical resonators via trapped ions 0, Cooper 
pair boxes 0, entanglement swapping [sl], and sudden 
switching of electrical interactions |9], have all been pro- 
posed. In the earliest proposal, [2], the entanglement 
provided a means for measuring the decoherence rate 
of coherent superpositions of nanomechanical resonator 
states. More recently, a scheme for entangling a super- 
conducting coplanar waveguide field with a nanomechani- 
cal resonator, via a Cooper pair box within the waveguide 
[10], was proposed. 

Here we propose a different scheme for entangling the 
nanomechanical resonator, based on the capacitive cou- 
pling of the resonator with the central conductor of the 
superconducting, coplanar waveguide, and which does 
not require any Cooper pair box (see Ref. [11] for a 
related proposal). The paper is organized as follows. 
In Sec. II we derive the Quantum Langevin equations 
(QLE) of the system and discuss when they can be lin- 
earized around the semiclassical steady state. In Section 
III we study the steady state of the system and quan- 
tify its entanglement by using the logarithmic negativity, 
while Section IV is for conclusions. 



cavity is driven at a frequency coo = cOc — ^o- In view of 
the equivalent circuit, the effective Hamiltonian for the 
coupled system is 
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where {x^Px) are the canonical position and momen- 
tum of the resonator, and Q) are the canonical co- 
ordinates for the cavity, representing respectively the 
flux through an equivalent inductor L and the charge 
on an equivalent capacitor C. The coherent driving 
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FIG. 1: Schematic representation of the capacitive coupling 
between a nanomechanical resonator and a superconducting 
coplanar microwave cavity, (a) Plan view of the device; the 
dark region is etched away and the resulting cantilever is met- 
alized to form one plate of a capacitor, (b) Equivalent circuit. 

of the cavity is given by the electric potential e{t) = 



-iujot 



). The function Co{x) de- 



scribes the capacitive coupling between the cavity and 
the resonator, as a function of the resonator displace- 
ment X. Expanding this around the equilibrium position 
of the resonator at d from the cavity and with capaci- 
tance Co, we have Co{x) = Co{l — x{t)/d). Expanding 
the capacitive energy as a Taylor series, we find to first 
order. 



II. QUANTUM LANGEVIN EQUATIONS AND 
THEIR LINEARIZATION 



The proposed scheme is shown in Figure [T] a nanome- 
chanical resonator is capacitively coupled to the central 
conductor of a superconducting, coplanar waveguide that 
forms a microwave cavity of resonant frequency oOc The 
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(2) 



where Cs = C + Co and /? = Co/Cs- 

We can now quantize the Hamiltonian, promoting the 
canonical coordinates to operators with [x,Px] = [Q,^] = 
ih. The quantum Hamiltonian, in terms of the raising 
and lowering operators for the cavity {a^ , a) and the res- 
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onator dimensionless canonical operators is 



where 
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and the couphng depends on 
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Typicahy, ujc ^ since uJc/^tt ^ 10 GHz [12], while 
cOm/^TT ~ 20 MHz [l3|. It is convenient to move into an 
interaction picture with respect to hw^a^a^ and neglect 
terms oscillating at ±2c(;o. The resulting Hamiltonian is 

Hi = hAoa'^a + ^^{f + p^) + hGoqa^a - ihE{a - a^). 

(7) 

The coupling term represents a low frequency modulation 
of the cavity resonance frequency. This will cause a phase 
modulation of the cavity field and write sidebands onto 
the cavity spectrum at multiples of cOm from coc- 

The resonator has a mechanical damping rate 7^ and 
the cavity bandwidth is 2k.. System dynamics also de- 
pend on the cavity input noise a'^^{t)^ where 



(a™'t(t)a»(t')) =n„5(t-t'), 



(8) 



with Ua = {exp{hu;c/kBT} — 1)~^, and also on the Brow- 
nian noise acting on the cavity ends ^(t), with correlation 
function llj 
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Clearly, ^(t) is not delta-correlated and does not de- 
scribe a Markovian process. However, quantum effects 
are achievable only when using resonators with a large 
mechanical quality factor {Qm = ^m/lm^ 1), and in 
this hmit ^(t) becomes delta-correlated |15j , 

mm + mm) n ^ 7™ i'^m + 1) 5{t - f), m 

where = {exp{hiJrn/kBT} — 1)~^, and we recover a 
Markovian process. Adding these inputs to the equations 
of motion that follow from ([7j), we obtain the nonlinear 
quantum Langevin equations (QLEs) 



q = UJrnP, 

p = -Umq-lmP^GQo)a^(^, 
a = — (/^ + iAo)a + zGoag + £^ ■ 
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Neglecting the noise and treating the deterministic equa- 
tions as classical, with a ^ a a complex field amplitude, 
we find the fixed points of the system by setting the left 
hand side of Eqs. (pT]) to zero. The fixed points are then 
given by 
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where the steady state photon number in the cavity is de- 
fined as ng = l^sp. Eq. ([T4j) is the same as the equation 
of state for optical bistability in a dispersive non linear 
medium [l^ and thus we expect for Aq > there will be 
multiple stable fixed points. 

The quantum dynamics of the full nonlinear system is 
difficult to analyze so we linearize around the semiclassi- 
cal fixed points. That is, we write a = 0^5 + (5a, q = qs-\-Sq 
and p = Ps -\- Sp. This decouples our system into a set 
of nonlinear algebraic equations for the steady-state val- 
ues and a set of QLEs for the fluctuation operators. The 
steady-state values are given by Eqs. (p!2|) and (p!3|) . and 
as = E/{k -\- iA); an implicit equation for a^, since the 
effective detuning A is given by A = Aq — Go|asP/^m- 
The QLEs for the fluctuations are 



Sq = Um^p, (15a) 
Sp = -u;mSq-^rn^p^Go{asSa^^a*5a) 

+(5a^(5a + ^, (15b) 
Sd = —{k ^ iA)Sa ^ iGo {as ^ Sa) Sq ^ V2K,a'^\lbc) 

Provided the cavity is driven intensely, \as \ ^ 1, we can 
safely neglect the terms Sa^Sa in Eq. (|15b|) and SaSq in 
Eq. (|15cp . and obtain the linearized QLEs 

Sq = cJmSp, (16a) 
Sp = -u;mSq-'ymSp^Goas{Sa^ ^Sa) ^£,,{16h) 

Sd = -{K^iA)Sa^iGoasSq^V2Ka'''^, (16c) 

where we have chosen the phase reference so that as can 
be taken as real. 



III. STEADY STATE OF THE SYSTEM AND 
ITS ENTANGLEMENT PROPERTIES 

In order to characterize the steady state of the sys- 
tem, it is convenient to rewrite Eqs. ([16]), defining G = 
Go^sV^, in terms of the field quadratures SX = {Sa + 
Sa'^)/V2 and SY = -i{Sa - 5a^)/V2, that is. 



Sq = Um^p, 

Sp -uJmSq - jmSp + GSX 
SX = -kSX^ASY^v^X' 
SY = -kSY - ASX ^GSq^ 



2kY' 



(17a) 
(17b) 
(17c) 
(17d) 
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where = {5a''' + ^a^^'tj/y^ and = -i{5a''' - 
5a'^''^) / ^/2. In matrix form, Eqs. (fT7|) can be written as 



u(t) = Au{t) + n{t), (18) 
{5q{t),5p{t),5X{t),5Y{t)), n^{t) = 



where 

(0,^(t), V^X^^(t), V^y^^(t)) and 
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Eq. (p!8|) has the solution 



(i5M(s)n(5), 



(20) 



where M{t) = exp(At). The stabihty conditions can be 
derived by applying the Routh-Hurwitz criterion [13], 

Si = 2^rnn I + {uJrn - + (^m + A) 



+7m [(7m + 2/^) {n^ + A^) + 2/^CJ^] } 
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+Acj^G^ (7^ 
cj^ (A^ + 



-2t^y > 0, 

G^A > 0. 



(21a) 
(21b) 



For driving on the blue sideband of the cavity (A = —ujrn) 
we have 

+27^^.3 + 7^c^^) } ^ / {^milm + 2/.)} , (22) 

while for driving on the red sideband (A = c^m), 

G < ^M^. (23) 

Since the noise terms in Eq. (p!8|) are zero-mean Gaus- 
sian and the dynamics are linear, the steady-state for 
the fluctuations is a two- mode Gaussian state, fully char- 
acterized by its symmetrically-ordered 4x4 correlation 
matrix. This has components Vij = {ui{oo)uj{oo) + 
Uj{oo)ui{oo)) /2. When the system is stable, using ([2Q|) . 
we get 



)ki, (24) 



k,i •'0 Jo 



where $(s - s')fc/ = {n{s)kn{s')i + n{s')in{s)k) /2 is the 
matrix of stationary noise correlation functions. Here 
*(s - s')ki = DkiS{s - s'), where D = Diag[0, 7m(2n;, + 
1), 2/t(na + 1/2), 2/t(na + 1/2)] and (gH) becomes 



V 



dsM{s)DM{sf 



(25) 



which, by Lyapunov's first theorem [l8|, is equivalent to 
AV + VA^ = -D. (26) 



Solving this equation, we can then quantify the entan- 
glement of the steady-state by means of the logarithmic 
negativity, Ej\f [H, [2^. This entanglement measure is 
particularly convenient because it is the only one which 
can always be explicitly computed and it is also additive 
[21I. In the continuous variable case we have 



Ej\f = max[0, — In 2r]] 

where r] = 2'^!^ [s(V) - ^{Vf - 4dety] 
S(y) expressed in terms of the 2x2 block matrix 



(27) 
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V 
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, with 



(28) 



as E(l/) = det 14 + det Va - 2 det Vcorr- 

The logarithmic negativity, assuming A = :^uOm^ de- 
pends on T, G, uoc^ ^, oom and 7^. We first consider the 
zero-temperature entanglement, such that our results are 
independent of T and uOc In all cases, entanglement in- 
creases with increasing coupling G; the limit on our en- 
tanglement being due to the limit on G specified by our 
stability conditions, ([22|) and ([23|) . and we shall set G just 
below this threshold. At zero temperature, the absolute 
magnitude of uOm is also insignificant, so we may hold it 
fixed (at uorn = 5 x lO^s"-^, say), leaving tz and 7^ as 
our remaining free parameters. It is implicitly assumed 
here that damping rates are controllable independent of 
resonant frequencies. The zero- temperature E_\f is shown 
in Figures 2(a) and 2(b) for driving on the blue and red 
sidebands, respectively. From this data, along with the 
stability conditions, we note that on the blue sideband, 
entanglement is maximized in a regime where cOrn ^ G, 
This is not the case for driving on the red sideband. We 
also observe that the logarithmic negativity plateaus in 
both cases as and 7^ increase. 

Now the temperature dependence of the entanglement 
follows from the Planck distributions specifying the noise 
input correlation functions, Eqs. ([37|) and (p!Q|) : hence the 
magnitudes of the resonant frequencies become signifi- 
cant in these calculations. Typical temperature depen- 
dence of the logarithmic negativity is shown in Figure 
3(a), decreasing from a positive value at zero tempera- 
ture to zero at a temperature we shall refer to as the 
critical temperature, Tc. Now Tc increases both with in- 
creasing ujc and uJrn'', henceforth, we shall consider these 
fixed, with coc = 10^^ s~^. The dependence of Tc on the 
damping is shown in Figures 3(b) and 3(c); the entangle- 
ment in the red sideband case appears more robust with 
respect to increases in temperature. 



A. Results in the rotating wave approximation 

The results in the regime cOm ^ G^n may be under- 
stood with the aid of a rotating wave approximation 
(RWA) calculation. We will refer to this as weak driv- 
ing as G depends on n^, the intracavity photon num- 
ber, which itself depends on the strength of the driving 
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FIG. 2: Zero-temperature logarithmic negativity for (a) driving on the blue sideband, (b) driving on the red sideband, n — 

K X 10-^S-\ 7m = 7m X IQ-^S-^ 




FIG. 3: (a) Typical temperature, T (mK) dependence of the logarithmic negativity (in this case, driving on the blue sideband 
with ^ = 7^ = 5 X lOS-^). (b) Tc (mK) as a function of damping for driving on the blue sideband, (c) Tc (mK)as a function 
of damping for driving on the red sideband. R — k x 10~^s~^, 7m — 7m x 10~^s~^ 



field, E. It is then useful to introduce the nanomechan- 
ical annihilation operator 5b = {Sq + i6p) / such that 
Eqs. ()16ap and (|16bp are equivalent to 

5b = -icUmSb-^ {6b - ((5at + ^a) + -^ (29) 

and the whole system is described in terms of annihi- 
lation and creation fluctuation operators by Eqs. ()16cp 
and ([29|) . We now move to a further interaction pic- 
ture by introducing the slowly-moving tilded operators 
db{t) = Sb{t)e-'^-^^ and da{t) = 5a{t)e-'^K They obey 



the QLEs 

si = _ll!L(sb-Sb^e^'^-')^^^-^, 
2 V J 

+V^a^^e^^^ (31) 
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The RWA allows us to ignore terms rotating at ±2cj^ 
and use A ^ Aq. Then, for driving on the blue sideband 
we have coc — ooo — oOm and 



6a 



2 



and for driving on the red sideband we have uOc 
and 



5h 



5a 



2 2 
-K5a + i-^^b + 



(32) 
(33) 

(34) 
(35) 



Note that a^^{t) = a*^(t)e*^^, possessing the same cor- 
relation function as a^^(t), and = ^(t)e^^-Vv^ 
which, in the limit of large c^m, acquires the correlation 
functions [22[ 

= ni,6{t-t'), (36) 
(6^"(t)6^"'t(t0) = [n^^l]5{t-t'). (37) 

From Eqs. ([32|) -(|33 |) we see that, for driving on the blue 
sideband, the cavity mode and nanomechanical mode 
play the role of the signal and the idler of a nondegener- 
ate parametric oscillator, characterized by an interaction 
term 5¥ 5a) + 5a5h. Therefore, it can generate entangle- 
ment. However, from Eqs. (|34|) - ([35]) . in the red sideband 
case the two modes are coupled by the beamsplitter-like 
interaction ^b'^^a + 8a) 8h^ which is not able to entangle 
modes starting from classical input states [23]. 

Now introduce tilded quadrature operators 5X = ((5a + 
5b))/^/2 and 5Y = {6a — 5d))/i^/2^ with corresponding 



= {a'' 



't)/y2, = 



input noise operators X* 

(a^^-a^^'t)/zv^, q^^ = (6^^ + 6^^'t)/^/2, andp^^ = (6^^- 
6^^'t)/iv^. We again obtain a system of the form (p!8|) . 
now with u^{t) = {6q{t),5p{t)),6X{t),SY{t), n^{t) = 
iV^q'^'it). V7^P^"(^), V^X*-(t), V2^Y'^{t)) and 



A = A^ = - 
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(38) 



where the upper (lower) sign corresponds to the blue 
(red) sideband case. For driving on the blue sideband, 
the stability condition of Eq. (|21ap simplifies in the RWA 
limit to 



G< ^2^ 



7^, 



(39) 



while the system is unconditionally stable for driving on 
the red sideband. For the symmetrically-ordered corre- 
lation matrix, we obtain an equation of the form 
which can be solved to give a matrix of the form 




\vt o" vtj 



V = v^ 

















(40) 



where 



^1^ 



rib 



1 2G^K [{na + 1/2) ± (nb + 1/ 2)] 

2 (7™ + 2k) (27™k t G2) 



{,41a) 



33 



1 , G V [(Ub + 1/2) ± (ng + 1/2)] 

2 (7™ + 2«)(27„ktG2) > 



2G-tmK [{fib + 1/2) ± (Ug + 1/2)] 

(7™ + 2k) (27„k t G2) 



(41c) 



Now deiV^^^ = tC^m)^^ which is non-negative in the 
red sideband case, a sufficient condition for the separa- 
bility of bipartite states [IJ]. Thus, in the red sideband 
case (RWA regime), the steady-state is not entangled. 

We can quantify the entanglement by proceeding along 
the lines of ([27j) and ([28|) . We may reproduce the entan- 
glement of Figures 2 and 3 for the blue sideband case, 
but we see no entanglement for the red sideband case. 
This is because the RWA regime puts us at a coupling 
far below the instability threshold. When the blue side- 
band steady-state correlation matrix is symmetric (that 
is, 2k = ^rn and Ua = Ub = ^) we find 



Ej^ — max[0. In 



1^G/2k 
l + 2n ^ 



(42) 



This and the stability condition of Eq. (|39|) imply that 
entanglement vanishes when 2n > 1, and that the loga- 
rithmic negativity is bounded above as Ej^f < ln2. Com- 
parison with Figure 2(a) shows that this is actually an 
upper bound in all cases. 

We shall now consider the experimental accessibility of 
the parameters described above. The coupling of Eq. (|6]) 
may be calculated by assuming (3 = 0.002, d = 100 nm, 

kg, 



6 X 10^° s-\ and m = 10" 



LJm ~ 10^ S"\ CJc 

giving Go ~ 18 s~^. For these parameters the equivalent 
capacitance is G ^ O.lpF and an equivalent inductance 
of L ^ 3nH. 

For driving on the red sideband and the largest damp- 
ing considered (7^ = k 

G < ^UjI^ -h k'^ ^ COm, 



= 10^ 
so G 



s ^), stability requires 
< 10^ s"^ Maximal 



coupling, before loss of stability, then corresponds to 
as = 10^, or a peak driving potential of 136 mV, which 
would be feasible. 

For driving on the blue sideband, the stability condi- 
tion in Eq. ([22|) . with k = 7^ and cOm ^ ^ stability 
requires G < \J2k ~ 10^s~^. Maximal coupling then 
corresponds to corresponding to olq — 10^, and a cor- 
responding maximum voltage of 38 /iV. It should also 
be noted that the very lowest damping rates depicted in 
Figures 2 and 3 would not be achievable, due to the finite 
quality factors of the resonator and cavity. 



IV. CONCLUSIONS 

We have shown a scheme able to entangle at the steady 
state a nanomechanical resonator with a microwave cav- 
ity mode of a driven superconducting coplanar waveg- 
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uide. The nanomechanical resonator is capacitively cou- 
pled with the central conductor of the waveguide and the 
steady state of the system, in an appropriate parameter 
regime, is entangled up to temperatures of tens of mil- 
liKelvin. We have explained how this can be achieved by 
presenting an approximate treatment based on a rotating 
wave approximation. 

Let us briefly discuss how to detect the steady state 
entanglement. From the above equations, and especially 
the correlation matrix of the steady state in the RWA 
limit, Eq. (|4Q|) . it is clear that the entanglement appears 
as a correlation between Sq{t) and SY{t), and also as a 
correlation between Sp{t) and SX{t). A measurement of 
entanglement thus requires that we measure these corre- 
lation functions. This is not an easy matter as it will re- 
quire highly efficient measurements of both the nanome- 
chanical resonator displacement and the field amplitudes 
in the microwave cavity. Methods based on single elec- 
tron transistors now enable a displacement measurement 
at close to the Heisenberg limit [13| . Unfortunately mea- 



surements of the weak voltages on the coplanar cavity 
are not yet quantum limited due to the need to amplify 
the signals prior to detection. This is not a fundamen- 
tal problem and a number of efforts are underway to 
do quantum limited heterodyne detection of the cavity 
fields. It thus seems likely that a direct measurement 
of the entanglement between a mesoscopic massive ob- 
ject and an electromagnetic field may be demonstrated 
using the approach of this paper. This would provide a 
path to entangling many nanomechanical resonators via 
a common microwave cavity field. 
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